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ABSTRACT 

We allow a more general (step-function) form of the primordial power spectrum 
than the usual featureless power-law Harrison-Zeldovich (with spectral index n = 1) 
power spectrum, and fit it to the latest Cosmic Microwave Background data sets. Al- 
though the best-fitting initial power spectrum can differ significantly from the power- 
law shape, and contains a dip at scales k ~ 0.003 hMpc~ , we find that fl m ~ 0.24, 
consistent with previous analyses that assume power-law initial fluctuations. We also 
explore the feasibility of the early releases of the 2dF and SDSS galaxy redshifts sur- 
veys to see these features, and we find that even if features exist in the primordial 
power spectrum, they are washed out by the window functions of the redshift surveys 
on scales k < 0.03 /iMpc -1 . 
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1 INTRODUCTION 



With the release of new Cosmic Microwave Back- 



ground (CMB) data from DASI ( Halvcrson ct al. 2001 \ 



ample, in supersymmetric double inflationary models, with 
two inflaton fields and two 'trigger' fields, the power spec- 
trum of density fluctuations is found to have a bump with 



BOOMERANG dNctterfield et al. 200 lD and MAXIMA flLcc 



et al. 2001), and the 2dF redshift survey (Percival et al 



superimpo sed oscillations on intermediate scales (Lesgour- 
gues 2000 ). Features in Pi n (k) can also be produced in dou- 



2001) nearing its completion, our ability to constrain cos- 
mological models has improved significantly. The new CMB 
data re moved the 'baryon crisis' caused by the unexpectedly 



low am ffitude of the second peak in the temperature power 



ble inflation (Polarski & Starobinsky 1992), and in one-field 
inflation with a feature in the inflaton potential (Starobin- 
sky 1992). The fluctuation spectrum may be se nsitive to 
physics at length scales below the Planck length (|Branden- 



spectrum, and the standard model of the universe now seems 
to be a flat Friedmann-Robertson- Walker universe, with 30% 
matter (~ 5% baryons + ~ 25% cold dark matter (CDM)), 
70% is dark energy, commonly parameterized by a cosmo- 
logical constant, and with the current value of the Hubble 
parameter Ho being around 70 km sec -1 Mpc -1 (e.g. Efs- 
tathiou et al. 2001; Wang, Tegmark & Zaldarriaga 2001). 

So me assumptions about the underlying coamological 
model are necessary in order to extract these parameters 



berger fc Martin 200 ij [Martin fc Brandenberger 200 1| ), and 
attempts have been made at extracting mass fluctuation 
spectra from models inspired by string theory (Khoury et 
al. 2001; Kempf 2001; Kempf & Niemeyer 2001; Easther, 
Greene & Shiu 2001). In Kempf & Niemeyer (2001) and Eas- 
ther et al. (2001) the primordial power spectrum was found 
to be of the Harrison-Zeldovich form, with a normalization 
depending on a short distance cutoff. More realis tic models 
will prob ably give rise to a fc-dependent imprint (Easther et 



al. 2001). It has been pointed out that models of this type 



from the data. One common assumption is that the initial 
power spectrum of the density fluctuations is a featureless 
power law Pi n (k) oc k n , where k is the comoving wavenum- 
ber. The spectral index n is found empirically to be close 
to the Harrison-Zeldovich value n = 1. This scale invariant 
primordial power spectrum is what typically comes out of 
models for inflation. However, there is no definite model for 
inflation, and some models predict features in Pi n (k). For ex- 



may run into problems in the form of an excessive creation 
of Planck energy particles at the present era (Starobinsky 
2001), but it should at any rate be clear that the theoreti- 
cal motivation for investigating more general forms for the 
initial power spectrum of density fluctuations is substantial. 

On the observational side, several claims of indications 
for features in Pi n (k) have been made (Broadhurst et al. 
1990; Griffiths, Silk & Zaroubi 2001; Atrio-Barandela et 
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al. 2001; Barriga et al. 2001; Hannestad, Hansen & Vil- 
lante 2001; Einasto et al. 1999; Gram ann fc Hiitsi 2001; 
Silberman et al. 2 001). Griffiths et al. (2001) and Hannes- 
tad et al. (2001) found that the CMB data favour a 
bump-like feature in the power spectrum at a scale fc ~ 
0.004 /iMpc -1 (h is the dimensionless Hubbl e para meter: 
H = lOO/i kmsec" 1 Mpc" 1 ). Barriga et al. fl200lh intro- 



duced a step-like feature in the range fc ~ 0.06-0.6 ftMpc -1 
and found that this spectral break gave a good fit to both 
the CMB data and the data from the APM survey (Maddox 
et al. 1990). Atrio-Barandela et al. ( 2001 ) investigated the 
temperature power spectrum in CDM models with a mat- 
ter power spectrum P m {k) at redshifts z ~ 10 3 of the form 
P m (Jfe) ~ fc" 1 ' 9 for fc > .05 ftMpc" 1 . This form was derived 
by Einasto et al. ( 1999 ) by analysing observed power spec- 
t ra of galaxies and clusters of galaxies. Gramann & Hiitsi 
(2001) studied the mass function of clusters of galaxies with 
this form of Pi„(k) and found that the predicted number 
density of clusters was smaller than the observed one. How- 
ever, these authors found that they could get a good fit to 
the mass function with a Pi„(k) having a dip-like feature 
at fc ~ 0.1 /iMpc" 1 , and that this Pin(fc) also was consis- 
tent with data from other cosmological probes like peculiar 
velocities and CMB. 

One of the reasons for considering alternatives to a 
scale- invariant Pi„(k) was that the CMB data before April 
2001 indicated that the amplitude of the second acoustic 
peak in the temperature power spectrum was low, result- 
ing in a baryon density Q, b h 2 ~ 0.03 outside the limits set 
by Big Bang Nucleosynthesis (BBN), which gives a 95 % 
confidence interval Q. b h 2 = 0.020 ± 0.002 (Buries, Nollett 
& Turner 2001). The new CMB data show a higher second 
peak, and the values for Q, b h 2 obtained with a power-law 
Pi„( k) a re now consistent with standard BBN ( Wang et al 
200lF" 



The motivation for the work presented in this paper is 
different: we want to check if the presently available CMB 
data allow for deviations from a scale-invariant Pi n (k), and 
if this is the case, whether these can be seen in the early re- 
leases of the 2dF and SDSS data. We will therefore consider 
more general shapes for Pi„(k) in our analysis. 

For simplicity, we will in our analysis vary only Q. m and 
Pin(k). We assume a flat Universe, and consider parameters 
like fi;, and h to be well constrained by other cosmological 
probes. The reader should note that these rather restrictive 
assumptions mean that the error bars we obtain for the es- 
timated quantities will be much smaller than obtained in an 
analys is with more free parameters, like that of Wang et al. 
(2001). Our aim in this paper is not to do the most general 
model fitting of the CMB+2dF, but to check two things: 
firstly, do the current CMB data allow features in Pi n (k) ? 
Secondly, can these features be seen in the early releases of 
the 2dF and SDSS data ? 

The outline of this paper is as follows: In Section 2 we 
introduce the models for Pi n (k) under consideration, and 
describe our procedure for fitting them to the CMB data. 
In Section 3 we present the results of this procedure, and 
in Section 4 we discuss possible constraints from the 2dF 
and SDSS galaxy redshift surveys. Finally, in Section 5 we 
summarize and discuss our results. 



2 THE PRIMORDIAL POWER SPECTRUM 
AND CMB ANISOTROPIES 

The power spectrum of fluctuations can be written as 

P(fc) = P m (fc)T 2 (fc), (1) 

where T(fc) is the transfer function (which modifies the ini- 
tial power spectrum during the radiation dominated era). 
To investigate more general forms for P% n {k), we let 



Pin(fc) = AkS(k), 



(2) 



where A is a constant and S(k) parameterizes the deviations 
from scale invariant initial fluctuations, and set up the mod- 
els as follows. We modified the p ublicly available CMBFAST 
code (Seljak & Zaldarriaga 1996) to include two alternatives 
for S(k): 

• a 'sawtooth'-shape, with 'teeth' equally spaced in ln(fc). 

• a set of 'top-hat' steps, equally spaced in ln(fc) and with 
amplitudes di,i = 1, . . . , N 

To be specific, we d efined the 'sawtooth' spectrum following 
Wang & Mathews fl200C|) as 



S(k) = 



01, 

kj-k 
ki-ki_ 



Oi-l 



+ 



fci-fci- 



fc fcl — fcmin 

ki-i <k < h (3) 



where 



(i _l)/(jV_l) 



,N 



(4) 



with fc min = 0.001 /iMpc" 1 , and fc max = 0.1 /iMpc \ Since 
the connection between the harmonic £ in the CMB power 
spectrum Ce and fc is roughly I ~ kd,A where for a flat 
universe the angular-diameter distance to the last scattering 
surface is well approximated by (Vittorio & Silk 1991): 

where c is the speed of light, the wave-number range 0.001 < 
fc < 0.1 /iMpc -1 corresponds for Q m = 0.3 approximately 
to 10 < I < 1000 which is nicely covered by the CMB data. 

Alternatively, we let Pi„(k) be defined by a set of 'top 
hat' steps, 



(6) 



{ai, fc < fci 
a», ki-i < k < ki 
ajv, fc > fcjv_i, 

with ki given by (^). The comoving wavenumber fc is mea- 
sured in units h Mpc -1 , where h is the dimensionless Hubble 
parameter. 

Both spectra are thus completely specified by N and the 
values of ai, . . . , ajv. In our calculations we chose N = 4. We 
also let the matter density Q. m be a free parameter, but the 
other parameters were kept fixed at Hp = 72 km s~ 1 Mpc~ 1 
flFreedman et al. 200l[), Q. b h 2 = 0.02 ([Buries et al. 200l[), 



N v — 3.04 (see e.g. Bowen et al. 2001), and we assumed a 
flat universe with no massive neutrinos and no reionization. 
We were therefore left with a five-dimensional parameter 
space to search for the best-fitting model. This was done by 
minimizing the \ 2 statistic given by 



X 



: (p) = ^[AT- 2 - A7f (p)](C-%[A7- 2 - AT/(p)], (7) 
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Figure 1. The best-fitting 'sawtooth' and 'top hat' models com- 
pared with the data. For comparison, we also plot a model where 
we have fixed n = 1, f2 m = 0.24, and the other parameters (S7;,, 
h etc.) at the values given in the text, and with COBE normal- 
ization, which is higher than most of the new data points. 



where p = (O m , ai, 02, 03, 04) is the parameter vector, A7^ 2 
are the measured CMB fluctuations, 



AT?(p) = |L £ + i)c,(p), 



(8) 



C e (p) is the COBE-normalized output from CMBFAST, 
To = 2.726 K, Wu is the window function, and Cy is the co- 
variance matrix of the observations. The observations, A%, 



Wu, and CVj are taken from (Wang et al. 2001). (The ab- 



sence of the usual 1 /^-factor in Eq. (g) is a result of their 
definition of the window function). We computed the \ 2 ° n 
a grid of 10 5 models, and located the region of parameter 
space containing the global minimum. In this region we per- 
formed a more accurate search for t he optimal param eters 
using the downhill simplex method (Press et al. 1992). 



3 RESULTS 

In Fig. [j] we show the best models for the two types of Pi n (k) 
we consider. For comparison we have also plotted the power 
spectrum for a model with Pi n (k) = Ak (i.e. S(k) = 1 and 
Q m — 0.24). The best-fitting 'top hat' and 'sawtooth' mod- 
els both have \ 2 ~ 32 for the 24 data points. The high x 2 
values are partly caused by the bandpowers centered at t — 2 
and £ = 50 in the compilation of Wang et al. Removing these 
points leads to lower values for the minimum x 2 , but has no 
significant effect on our estimated values for Q m and the pa- 
rameters of P in (k). Note that the CMB only constrains the 
quantity fi m /i 2 , but since we have fixed h in our analysis, 
we obtain an estimate for Q m directly. In Fig. ^|we show the 
S(k) which provide the best fit to the data, with ai scaled to 
1 for clarity. The best-fitting parameters, mean values and 




Figure 2. The best-fitting 'sawtooth' and 'top hat' S(k). Vertical 
lines indicate the error bars on the estimated parameters. The k 
values spanned in this figure corresponds roughly to values of I 
in the range 1-10 4 . 



confidence intervals can be found in Table 1. The confidence 
intervals for the parameters were obtained by the standard 
approach of constructing marginalized likelihoods for each 
parameter by integrating out the other parameters from the 
likelihood C oc exp(— x 2 /2). The tight confidence intervals, 
in particular on fi m , reflect the restrictive assumptions we 
have made about other cosmological parameters. 

We see from Fig. |^ that both the best-fitting S(k) have 
a dip at k ~ 0.003 /iMpc -1 . In fact, for the 'top hat' Pi n (k) 
the data favour ai ~ a 2 and 03 ~ 04. This means that there 
are effectively only two parameters in this spectrum: the 
position of the dip and the relative sizes of the amplitudes. 
We have checked this by repeating the analysis with more 
amplitudes in the spectrum, and found that the data in this 
case still favour a Pi n (k) with a dip at k ~ 0.003 /iMpc -1 . 

Thus, in this case we can reduce the parameter space to 
three dimensions: f2 m , the ratio 1Z between the two ampli- 
tudes defining Pi n (k), and the position in fc-space k* of the 
dip. The best-fitting values and confidence limits are given 
in Table 2. For completeness, we also made a calculation for 
the 'sawtooth' spectrum with this reduced set of parame- 
ters (i.e. keeping just one 'tooth' in the spectrum, and us- 
ing its amplitude and position in fc-space as free parameters, 
see the tables) . The marginalized likelihood distributions for 
Q. m , Ti, and fc* are shown in Fig. |^. We see that Q, m is well 
constrained to a narrow range around 0.24, and the size of 
the dip 1Z similarly constrained to be around 1.6, consistent 
with the results for four steps. The scale fc* at which the 
break occurs has a broader distribution. In the calculation 
with four steps, this scale was at ~ 4.6 ■ 10 -3 /iMpc -1 , but 
from Fig. | and Table 2 we see that when we allow fc* to 
vary, we can only constrain it to be in the range ~ 0.001- 
0.005 /iMpc -1 . 



4 COULD 2dF AND SDSS DETECT 

FEATURES IN THE PRIMORDIAL POWER 
SPECTRUM ? 

Combining data from different cosmological probes can in 
many cases lead to tighter constrain ts on the cosmological 
param eters, see e.g. Efstathiou et al. ( 2001 ) and Wang et al. 
(|200lj ). Since we saw in the previous section that the CMB 
does not rule out deviations from a scale- invariant Pi„(fc), 
it is interesting to see if we can obtain further constraints 
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Figure 3. Marginalized likelihood distributions for the parame- 
ters in the 'top hat' model. 



Table 1. Best- fitting values and 95 % confidence intervals for the 
parameters of the two models. 





'Top hat' 
Best fit 


95% c.i. 


'Sawtooth' 
Best fit 


95% c.i. 




0.24 


24+ ' 02 


0.29 


29+ 02 


ai 


0.56 


59+ 012 

u - oa -o.io 


1.0 


1 00+ 02 


ai 


0.55 


59+ ' 11 


0.57 


57+ 013 


a-3 


0.38 


40+ ' 04 


0.74 


70+ ' 07 
u ' '-0.12 


a 4 


0.37 




0.47 


0.38±°;S 



from the matter power spectrum as estimated from the 2dF 
and SDSS galaxy redshift surveys. 

Assuming a simple scale-independent biasing model 
with a bias parameter b, the galaxy power spectrum (lin- 
ear theory) is predicted to be 



P S (k) 



b 2 P 

U JTrr 



(ft) = b' 2 AkS(k)T 2 (k) 



(9) 



where T(k) is the transfer function. The recent analysis of 
2dF+CMB data by Lahav et al. (2001) found b = 1.0±0.1 on 
comoving scales of 0.02 < k < 0.15 ftMpc^ 1 . Using CMB- 
FAST, we computed P(k) for S(k) = 1 and the best-fitting 
'top hat' and 'sawtooth' S(k), with the results shown in Fig. 
^| However, to compare with the galaxy power spectrum 
from 2dF (Percival et al. 2001), we must convolve P g (k) 
with the 2dF window function: 



/(fc) tx 



|k-qixiwgV<z, 



(10) 



where { | | 2 ) is the spherical average of the the 2dF window 
function, approximately given by 



Table 2. 95 % confidence intervals for the parameters of the 
reduced models with two bins. As in Table 1, the central values 
are the mean values, k* is given in units of 10 — 3 /iMpc -1 . 





'Top hat' 




'Sawtooth' 






Best fit 


95% c.i. 


Best fit 


95% c.i. 


TZ 

k* 


0.24 
1.55 
2.62 


24+ - 02 

1 63+ ' 21 
i.oo_ . 2g 

2 9+ 1 ' 6 
z - y -2.0 


0.24 
1.74 
5.0 


n 9 o+0.03 
u - zo -0.03 
i 77 +0.26 
'-0.34 
4 7+3-0 
^■'-4.1 
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Figure 4. Theoretical power spectra from CMBFAST. The solid 
line corresponds to S(k) = 1 (H m = 0.24, h = 0.72), the dashed 
line to the best-fitting 'top hat' S(k) model, and the dash-dotted 
line to the best-fitting 'sawtooth' S(k) model. 



(\w k \ 2 ) = 



1 + (fe/a) 2 + (k/b) 4 



(11) 



with a = 0.00342, and b = 0.00983. As shown in Appendix 
A, the convolution integral can be rewritten as 



(ft) cx / K(k,k')P g (k')dk', 



(12) 



where K(k, k') is given by an integral over the 2dF win- 
dow function ( |ri"| ) that can be evaluated analytically. As 
discussed in Appendix A, for k ~ 0.1 /iMpc -1 , K(k,k') is 
almost a delta function 8(k' — k), so in this region P CO nv(k) = 
P(k). For low k, K(k,k') is a broad distribution, and the 
main contribution to the convolution integral comes from 
values of k! larger than ~ 0.01 ftMpc" 1 , so that P C onv(fc) is 
nearly independent of k. 

As a result, the features in P(k) introduced by S(k) are 
washed out by the convolution, as can be seen from Fig. ^. 
The results show that the present data from the 2dF sur- 
vey cannot give us information about the power spectrum 
at wavenumbers smaller than about 0.03 /iMpc -1 , every- 
thing on larger scales is washed out. There is no relation 
between the wiggles visible in 2dF power spectrum, which is 
the result of observing a single realization of the true power 
spectrum convolved with the window function of the survey, 
and the features we introduced in Pi„(k). The wiggles in 
the 2dF power spectrum may be signatures of baryon os- 
cillations, but may also be a result of correlated noise. The 
scale of the observed wiggles is Aft = 0.03 Mpc -1 , which is 
the same scale over which power is correlated, so the wig- 
gles could well be the result of correlated noise. This ten- 
tative conclusion is supported by the recent analysis of the 
published sampled carried out by Tegmark, Hamilton & Xu 
(2001). 

We also compared our calculated P g (k) with the recent 
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Figure 5. Theoretical matter power spectra convolved with the 
2dF window function and compared with the 2dF data. The solid 
lines corresponds to S(k) = 1 (f2 m = 0.24, h = 0.72), the dashed 
line to the 'top hat' S(k), the dash-dotted line to the 'sawtooth' 
S(k), and the vertical bars are the data points from Percival et 
al. (2001). 



estimate of the power spectrum from the S loan Digital Sky 
Survey (SDSS), taken from Dodelson et al. (2001), where the 



SDSS data were analysed to obtain the deconvolved three- 
dimensional power spectrum. The SDSS results are given 
as a set of bandpowers in fc-space, and thus the calculated 
power spectrum P g (k) must be transformed in a way analo- 
gous to Eq. (^); in the ith bin, the bandpower P(fa) is given 
by 



E 



W S dss(», 3)P B {kj). 



(13) 



dow functions Wsdss can be found in Dodelson et al. (2001). 
Our results are shown in Fig. |(| Only results for the magni- 
tude bin r* = 21 — 22 are shown, but the same conclusion is 
true for the other three magnitude bins given in Dodelson 
et al. (2001). We see that the same conclusion applies to the 
SDSS spectrum as to the one from 2dF: at present there is 
no information on the scale where we find features in P;„(fc). 
However, this conclusion only applies to the presently avail- 
able data, as both 2dF and SDSS will, once completed, give 
information about larger scales than those probed by the 
data used in this paper. 



5 CONCLUSIONS 

We have analysed the new and updated CMB data, relaxing 
the usual assumption of power-law initial fluctuations. We 
found that the value of Q. m we could extract was consistent 



with other analyses, e.g. Sl„ 



0.24 



+0.02 



for the 'top hat' 



Pi„(k) with varying position of the break. We should point 



Figure 6. Comparison of the 'top hat' P g (k) with the SDSS 
power spectrum in the magnitude bin r* = 21 — 22. Note that, in 
contrast with Fig. 5, it is the deconvolved power spectra that are 
shown in this figure. 



out that the small error bar on Q m found in our analysis 
is mainly a result of the restrictive priors we put on other 
quantities. For example, allowing the Hubble constant to 
vary would have led to a significant increase in the uncer- 
tainty in Q m since the CMB power spectrum depends on 
S7 m only through the physical matter density u m = Q m h 2 . 

We find that the present CMB data allow the initial 
power spectrum of the density fluctuations, Pi n (k) to have 
quite significant features, in particular a dip at a comoving 
wavenumber k ~ 0.003 /iMpc -1 , which corresponds to t ~ 
40. This dip can be understood as follows: Increasing Q. m h 2 
decreases the amplitudes of the peaks in the CMB power 
spectrum, and also shifts their positions to lower multipoles 
£. The position of the first peak is well determined by the 
data, and the models pay a higher price in terms of the % 2 f° r 
not fitting it than they do for not fitting the amplitude. The 
position of the first peak in the compilation of CMB data we 
have used is well fitted by Q m h 2 = 0.12 (Wang et al. 2001), 
which for h = 0.72 gives n m = 0.24. Fig. |l| shows that the 
CMB power spectrum of a model with this value of Q m h 2 , 
but a Harrison-Zeldovich P;„ (k) , lies consistently above the 
data points for all but the lowest values of £. The freedom 
in the 'sawtooth' and 'top hat' Pi n (k) allows the amplitude 
to be fitted by reducing Pi„(k) at comoving wavenumbers 
above k ~ 0.003 hMpcr 1 . 

From the preceding discussion it also follows that the 
parameters of the 'sawtooth' and 'top hat' Pi n (k) are sen- 
sitive to calibration errors in the data points, and also that 
the choice of normalization for the models plays a role. We 
have normalized all our CMB power spectra to COBE, using 
the prescription given by Bunn & White (1997). Recent work 
by Lahav et al. (2001), Reiprich & Boehringer (2001), Seljak 
(2001), and Viana, Nichol & Liddle (2001), found that cr 8 , 
the rms mass fluctuations in 8 ft -1 Mpc spheres, is ~ 20% 
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lower than predicted by the COBE normalization. If we had 
used these results to normalize our CMB power spectra, the 
size of the dip we found in S(k) would have been reduced. 
However, our goal here was to test the maximal features 
that can still be consistent with the CMB data. Therefore 
the analysis with the constraint of COBE normalization can 
be viewed as an upper limit on the amplitude of the features. 

We also compared the early-release 2dF and SDSS 
galaxy power spectra to the theoretical predictions from our 
best-fitting Pi n (k), and we found that the observed power 
spectra are not sensitive to the features in Pi„(k) at comov- 
ing scales k < 0.03 ftMpc -1 . There is no relation between 
the features we find in Pi„{k) and the wiggles observed in 
the 2dF power spectrum. 

Furthermore, we point out again that we have fixed all 
parameters except Q m and Pi n {k) in our analysis. Combin- 
ing extra degrees of freedom in Pi n (k) with a full-scale anal- 
ysis of the CMB data, would lead to larger error bars on the 
cosmological parameters. The more accurate measurements 
of the CMB fluctuations expected in the near future, anal- 
ysed jointly with data sets from other cosmological probes, 
will hopefully allow us to put tighter constraints both on 
the primordial fluctuations and the parameters defining the 
geometry of the Universe. 



where 9 is the angle between k and q. (Here and in the 
following we ignore the normalization constant. In practical 
calculations it is taken care of by dividing by / (\W q \ 2 )d 3 q.) 
On substituting k' 2 — k 2 + q 2 — 2kqcos8, d(cos(9) = 
—k'dk'/kq, the integral becomes 



Pconv(k) OC — 

k 



dqq(\W q \ 2 ) 



k + q 



dk'k'P g (k') 



(A4) 



k — q\ 



With the simple approximation Eq. ( A2) to the 2dF window 
function, the convolution integral can be simplified further 
by changing the order of integration in Eq. (A4). We in- 
tegrate over the region in the (k' , g)-plane bounded by the 
lines k' — k + q, k' — k — q, q < k, and k' = q — k, q > k, so 
changing the order of integration is easy, and the result is 



dk'k'Pg(k') 



k' + k 



dqq(\W q 



K(k, k')P g (k')dk' 



where 
K(k,k') 



2nk' 
k 



k' + k 



dqq{\W q \ 2 ) 



(A5) 



(A6) 
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APPENDIX A: THE CONVOLUTION 
INTEGRAL 

The convolution of the power spectrum P(k) with the win- 
dow function is given by 



Pconv(k)cx / P 9 (k-q)|H/ k (q)| 2 d 3 g. 



(Al) 



Assuming P is isotropic, and that only the spherical average 
of the final power spectrum is of interest, one can replace 
the window function with its spherical average, which for 
2dF can be approximated by 



(\w k \ 2 ) 



1 



1 + ak 2 + /3k 4 



(A2) 



where a = 8.55 -10 4 , (3 = 1.071 -10 s . Choosing k as the z-axis 
in the integration, the convolution integral can be written 

as 

Pconv(fe) OC / P(|k-q|)(|W/„| 2 }d 3 9 



= 2tt / dqq 2 (\W q \ 2 ) 
Jo 
f+l 



X / d(cos 9)Pg ( yjk 2 + q 2 - 2kq cos (?) , (A3) 



With the 2dF window function, one can obtain an analytical 
expression for K(k, k'): 



K(k,k') = 



2-Kk! 



k' + k 

k'-k] 
(k' + k) 2 



dq 



1 + aq 2 + /V 
dx 



71-fc' 

k J (k ,_ k)2 1 + ax + flx 2 

= -^-?7{arctanh[£(fc' — k) 2 + A] 
k 

— arctanh[£(fc' + k) 2 + A]}, (A7) 

where <q = 1.2055 ■ 10 -5 , £ = 2.5822 ■ 10 3 , and A = 1.0307. 

Numerical plots of K(k, k') for various values of k shows 
that it can be roughly approximated by a Gaussian, 



K(k, k) oc exp 



(k' - fi k y 
1°l 



(A8) 



and that for k > 0.1 ftMpc" 1 , a k <C k, and fik ~ k, so that 
the Gaussian approaches 8(k' — k). In this regime of k, we 
will therefore have 



(k) cc / S(k' - k)P g (k')dk' = Pg(k). 



(A9) 



For k <gi 0.1 ftMpc -1 , the Gaussian is very broad, &u 3> k 
and crfe 2> /ifc. To illustrate what happens in this regime of 
k, we take P g (k) oc k for k < kb and P g (k) oc k~ 2 for k > kf,. 
Then 

Aonv(fc) OC f b exp (- ( ^ 7 ) /,',//,' 



+ 



exp 



(k'-» k ) 2 \ dti_ 
2o\ J k' 2 ■ 



(A10) 



Both integrals can be evaluated analytically in terms of the 
error function erf(s). Numerically, it turns out that it is 
a reasonable approximation in this regime to replace the 
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Gaussian by its peak value, which gives 

-fconv 

(k) oc 



Wang X, Tegmark M., Zaldarriaga M., istro-ph/0105091 



k'dk'+f 

Jk h 



Wang Y., Mathews G.,2000, astro-ph/0011351 



2 + k b 



1 

kb ' 



(All) 



Note that this is independent of k, which agrees with the 
numerical results for the convolved P{k). 



REFERENCES 



Atrio-Barandela F., Einasto J., Miiller V., Miicket J. P., Starobin- 

sky A. A., 2001, ApJ, 559, 1 
Barriga J., Gaztanaga E., Santos M., Sarkar S., 2001, MNRAS, 

324,977 

Ro wcn R... Hansen S H., Melchiorri A., Silk J., Trotta R., 2001, 



-ph/01 10636 



Koo D. C, Szalay A. S., 1990, Nat, 



6 

2001, 



Phys. Rev. D, 63, 



Shiu G., 2001, Phys. 



Brandenberger R. H., Martin J., 2001, Mod. Phys. Lett. A, 16 
999 

Broadhurst, T. J. Ellis R. S 
343, 726 

Bunn E. F., White, M., 1997, ApJ, 480 
Buries S., Nollett K. M., Turner M. S.. 

063512 

Dodelson S. et al., 2001, ^stro-ph/0107421 
Easther R., Greene B. R., Kinney W. H 

Rev. D, 64, 103502 
Efstathiou G. P. et al., 2001, a.stro-ph/0109152 
Einasto J. et al., 1999, ApJ, 519, 441 
Freedman W. L. et al., 2001, ApJ, 553, 47 
Gramann M., Hiitsi G., 2001, MNRAS, 327, 538 
Griffiths L. M., Silk J., Zaroubi S.. 2001 MNRAS. 324, 712 
Halverson N. W. et al., |astro-ph/0104489| 
Hannestad S., Hansen S. H., Villante F. L 

Phys., 16, 137 
Kempf A., 2001, Phys. Rev. D, 63, 083514 
Kempf A., Niemeyer J. C, 2001, Phys. Rev. D, 64, 103501 
Khourv J.. Ov rut B. A., Steinhardt P. J 



2001, Astropart. 



th/0 103231 



Turok N. 



hcp- 



Lahav O. et al., 2001, ^stro-ph/0112162 
Lee A.T et al., 2001, ApJ, 561, LI 
Lesgourges J., 2000, Nucl. Phys. B, 582, 593 

Lidsey J. E., Liddle A. R., Kolb E. W., Copeland E. J., Barreiro 

T., Abney M., Rev. Mod. Phys., 69, 373 
Maddox S. J., Efstathiou G. P., Sutherland W. J., Loveday J., 

1990, MNRAS, 242, 43P 
Martin J., Brandenber ger R.. H.. 2001. P hvs. Rev. D, 63, 123501 
Netterfield C.B. et al.. astro-Dh/0104460 
Percival W. J. et al., a.stro-ph/0105252 

Polarski D., Starobinsky A. A., 1992, Nucl. Phys. B, 385, 623 

Press W. W. et al., 1992, Numerical Recipes in Fortran, 2. ed. 
Cambridge University Press 

Reiprich T. H., Bochringcr H.. 200 1. 

Seljak U., 2001, |astro-ph/011136; 

Seljak U., Zaldarriaga M., 1996, ApJ, 554, 47 

Silberman L., Dekel A., Eldar A., Zehavi I., 2001, ApJ, 557, 102 

Souradeep T., Bond J. R., Knox L., Efstathiou G. P., Turner M. 
S., 1998, in Roszkowski L., ed, Proc. COSMO-97. World Sci- 
entific Proceedings of COSMO-97, Ambleside, England, Ed. 
L. Roszkowski, World Scientific 

Starobinsky A. A., 1992, JETP Lett., 55, 489 



astro-ph/0111285 



3. 371 



Starobinsky A. A., 2001, JETP Lett 

Tegmark M., Hamilton A. J. S., 2001, |astro-ph/0111575 
Vittorio N., Silk J., 1991, ApJ, 385, L9 

Viana P. T. P, Nichol R. C, Liddle A. R., 2001, |astro-ph/0111394 



© 0000 RAS, MNRAS 000, 000-000 



